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SUMMARY 


A method is outlined for calculating the expected number of maxima 
or m-in-imo of a random process with non-Gaussian frequency distribution 
from the statistical moments of the process and its first two derivatives. 
This method is based on an estimate of the joint frequency function of 
the process arid its first two derivatives given by means of a generalized 
form of Edgeworth’s series; the procedure thus consists essentially in 
applying a correction to the results for a Gaussian process. The func- 
tions required in this procedure are calculated for the first two correc- 
tion terms; therefore, the effects of skewness and kurtosis can be cal- 
culated, provided the required moments are known. Expressions are given 
for these moments in terms of multiple correlation functions and multi- 
spectra, and the relations between these functions for a random output 
of a linear system and those for the random input are indicated. 


INTRODUCTION 


Many physical processes of interest in aeronautics and allied fields 
are dete rmina te only in a statistical sense. Such processes are referred 
to as stochastic or random processes . If the statistical characteristics 
of such a process are invariant in time, it is referred to as a stationary 
random process. Hie basic problem in connection with these processes Is 
usually either to predict the output of a dynamic system which is subjected 
to a random inp ut (so that the output is also generally random in nature) 
from the statistical characteristics of the input and the dynamic charac- 
teristics of the system, or to estimate certain statistical characteris- 
tics of a given process from others. (See refs. 1 to 6 for discussions 
of several problems in communications theory and aeronautics from the 
point of view of random-process theory.) 

One statistical characteristic which is frequently of interest is 
the number of max ima or m-Tnimp. expected in a given time; that is, the 
number of positive or negative peaks of the process within a certain 
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range or exceeding a certain level that can be expected in that time. 

The expected life of- an airplane, for instance, depends on the expected 
number of times in a given period of time that its ultimate load is 
likely to be exceeded. (See refs. 4 and 5») Similarly, the fatigue 
life of a structure can in some cases be related to the number of maxima 
per unit time and their frequency distribution. (See ref. 6, for 
instance. ) 

For a stationary Gaussian process - that is, for a stationary random 
process in which the stochastic variable and its derivative are jointly 
normally distributed - Rice (ref. 1) has given a simple expression for the 
expected number of maxima in terms of the second moments of the process 
and its first two derivatives. These moments can be obtained from the 
correlation function or power spectrum of the process. In turn, if the 
process represents the output of a linear system, the spectrum can be 
related very simply to the correlation function or ..spectrum of the input. 

In the present paper a similar expression f or Jbhe expected number 
of maxima is obtained for a stationary process with a joint frequency 
distribution of the process and its first two derivatives; this distri- 
bution differs slightly from the normal. The approach used herein con- 
sists in expressing the joint frequency distribution of the process and 
its derivatives in terms of its second and higher moments by means of a 
multivariate form of Edgeworth's series, so that the desired expected 
number of maxima or minima can then be expressed in terms of these moments. 
Again, these moments can be expressed in terms of correlation functions 
and spectra, and the correlation functions and spectra of an output can 
be related readily to those of the input. However, in this case more than 
the ordinary (double) correlation function or power spectrum Is required, 
because the nth moments depend on the n-tuple correlation function or 
the corresponding spectrum. Hence, some of the multiple correlation 
functions or multispectra of the input must be known If the number of 
maxima of an output process with non-Gauss ian frequency distribution is 
to be predicted by this method. 

I nasmuch as the terms of Edgeworth's series represent, essentially, 
corrections to a normal distribution, the approach _outlined herein also 
furnishes, essentially, a correction to the results obtained for a 
Gaussian process. Explicit expressions are given herein for the func- 
tions required in the first two correction terms, which involve the third 
and fourth moments. No such expressions are given for higher correction 
terms, because the effort entailed in obtaining the required multispectra 
soon becomes very large . The procedure given herein furnishes an esti- 
mate of the effects of skewness and kurtosis on the results of interest 
and is, therefore, best suited to distributions which differ by relatively 
little from the normal one, that is, primarily in the third and fourth 
moments but to a lesser extent in the higher moments. 
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Consideration will be confined to a process with zero mean. This 
trivial restriction implies that, for a process with nonzero mean, the 
results given here apply directly to the process which is the difference 
between the actual process and its mean; they can, however, be modified 
to apply to the actual process in a straightforward fashion. 


SYMBOLS 


D n 


operator 



f(x 0 ,Xi,x 2 ) 

f G (x 0 ,Xi,x 2 ) 

F(t) 

h(t) 

H(co) 

I mnp ( x 0) 

J mnp ( x o) 

M ij 

n ( x o) 

N W 

t 

x(t) 

x(t),X(t) 


trivariate frequency-distribution function for a random 
process and its first two derivatives 

trivariate normal distribution (with zero means) 
random input process 

Indie ial-response function for linear dynamic system 

frequency-response function for linear dynamic system 

function used in expression for 
equation ( 17 ) 

function used In expression for 
equation (21) 

reciprocal second moments (elements of Inverse of 
matrix of second moments) 

expected number per unit time of maxima with inten- 
sities in band of unit width centered on xq 

expected number per unit time of maxima above xq 
time 

given random process, which may be output of a linear 
system subject to random input 

first and second time derivatives of x(t) 


n(x 0 j , defined in 
N(xq) , defined in 
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x 0> x l' x 2 


a 


random variables corresponding to x, x, and X, 
respectively 


parameter , 



a mnp 




joint statistical moments of random process and its 
first two derivatives, defined in equation (6a) 
or (6b) 

joint statistical moments of random process and its 
first two derivatives, defined in equation ( 5 ) 


A 


determinant of matrix of second moments 


i 

*x 

*xx 

>|r 

xxx 

<Px 

<Pxx 

^xxx 


1 xn 

dimensionless random variable, 

(double) correlation function for x(t) 

triple correlation function for x(t) 

quadruple correlation function for x(t) 

power spectrum for x(t), Fourier transform of i|r x 

double power spectrum for x(t), Fourier transform 
of *xx 

triple power spectrum for x(t), Fourier transform 

of W 


T time displacement, argument of ty x 

cd circular frequency 


ANALYSIS 


Basic Relations 

For a given stationary random process x(t), the number of maxima 
that are expected to exc eed the level Xq per unit time will be 
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designated by N^Xq^ , and the dumber of maxima in an intensity band of 
unit width centered on Xq that are likely to occur per unit time will 
be designated by n^XQ^ , so that 


n 



d X q 


( 1 ) 


In reference 1 the follow in g expression is given for 11 in 

terms of the joint frequency distribution f (xq,x-j_,x 2 ) of the process 
and its first two derivatives: 


n 



x 2| f ( x 0>°> x 2) dx 2 


( 2 ) 


where f (Xq,x-^,x 2 ^ is defined by the fact that f ^xq,x-^,x 2 ^ dxQdx-^dxg 
represents the probability that at time t: 

Xq < x(t) < Xq + dXQ 

X 1 = x(t) < x^ + dx-j_ 
x 2 S 5 f(t) < X2 + dxg 

This function is invariant with t by virtue of the assumed stationarity 
of the process. 

For a Gaussian process this result can be expressed in an especially 
simple form. For such a process the frequency function is 


f G (x 0 ,x 1 ,x 2 ) = i e -i( M 0O x O 2 -^ll x l 2 ^22 x 2 2+ ^02 x O x 2) (3) 

^ (2*) 3 / 2 \/A 


where A is the determinant of the second moments of the process, and 
the coefficients are the elements of a matrix reciprocal to the 

matrix of the second moments and will be discussed further in a later 
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section. Substitution of this expression for f (x Q ,x 1 ,x 2 ) in equa- 
tion (2) yields an expression for n(x Q ) in terms of the moments; hence, 
integration over x Q yields an expression for N(xq) , and thus also for 

total number of maxima per unit time, which is equal to N (-«>). The 
expressions for n^XQ^ and N(-<») are given in reference 1. 

For a non-Gaussian process a similar expression, for-the jo int 
frequency distribution in terms of the moments can be obtained from a 
multivariate version of Edgeworth's series. This series is derived (see 
ref. 7, for instance) on the assumption that the given process represents 
the sum of a large number M of statistically independent random vari- 
ables. Then, by expanding the characteristic function for the process 
in a series, several asymptotic expressions can be obtained (depending 
on the manner in which the terms are collected) for the given non- 
Gaussian distribution in terms of its moments, the Gaussian distribution, 
and its derivatives. As M tends toward infinity, all terms of the 
series except the one which represents the Gaussian ; part of the distri- 
bution tend to zero. In Edgeworth's series, terms are grouped according 
to powers of M, so that- each group can be expected to contain the terms 
representing a given extent of deviation from the Gaussian distribution. 

This derivation can readily be extended to multivariate distributions 
by using the concepts of random-vector theory, such as the multivariate 
form of the characteristic- function and of the Gaussian distribution. 

(See ref. 7, for instance.) The results can be expressed for the case 
of interest as _ 



li J 

fa» k S3 

[ J 

a ijkZ _ 5a ij a kZ 

L 1 

5! Bx.j_<jXjBx k 

1 

1 



+ 


lOa^V™ 



6 ! 


Bx^Bxj Bx^Bx 1 Bx m Bx n j 


• f G(*b> x l' x 2) ( 1 * * 4 ) 


1 ikl 

where a 0 *** represents the moments of the given process and the 

indices i, j, k, . . . may have the value 0, 1, or 2. In this expres- 

sion the summation convention is used, so that any index repeated in a 
product implies a summation over that index. Terms associated (in the 
derivation) with a given power of M are grouped within braces. Only 
the first two terms beyond the Gaussian part are listed here; from the 
Edgeworth's series given in reference J, one additional term can be 
deduced by analogy. 7 
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The Statistical Moments of the Process 
As used in equation (4), the moments are defined by 


ijk. .. = ~ 


( 5 ) 


where the bar designates a time average, and the symbols D , D , . . . 
designate the ith, jth, . . . derivative with respect to time. The 
reason for this definition is that it permits the application of the 
summation convention and thus greatly simplifies the writing of equa- 
tion ( 4 ), and the reason for the superscript notation is that subscripts 
will be used for the moments defined in the manner which is more descrip- 
tive and convenient for the purpose at hand, because it assigns only one 
set of indices to any moment, namely 


ctjjmp = x m (t)* n (t)X p (t) 


( 6 a) 


or 


/ 00 p CO p 00 

J x 0 m xi n x 2 P f(x 0 ,x 1 ,x 2 )dx 0 dx-L dx 2 


( 6 b) 


The two sets of moments defined by equations ( 5 ) and ( 6 ) can be 
identified with each other in the following manner: For an rth moment 

the number of superscripts i, j, k, . . . is r, and m+n + p= r; 
m, n, and p are, respectively, the numbers of 0 's, l's, and 2 ‘s among 
i, j, k, . . . . Thus, for instance. 


2 . ..3 7.' 001222 

yr x yr = 0213 mxxxxxx = a 

= a 222100 _ a 120220 

and so on. Obviously, the superscripts of a J can be permuted 

in any manner, but any change in the subscripts of cx^np changes the 
moment referred to . 

The first moments (mean values) of the process ccq_qo = x ’ ^lO = 
and = x are zero; x is zero by stipulation and the others are 

zero of necessity, inasmuch as the process is stationary. Of the nine sec- 
ond moments, four are zero and two are equal to each other as a result of 
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the assumed stationarity of the process; therefore, the following moments 
are left : 


“200 “ x2 
OQ 20 = x 2 


a 002 


Jf 2 


aid = x x 


As pointed out in reference 7 , for instance, they can all be expressed 
in terms of the (double) correlation function for -x(t), which is 
defined by — 


* X (T) = x(t) x(t+T ) 


( 7 ) 


or in terms of the power spectrum, which is defined by 





e" icOT 1 l; x (T)dT 


(8) 


The required relations are given in reference 7 and^are repeated for the 
sake of completeness in table 1. 


The second moments can be arranged in a matrix' in the following 
manner : 



a 200 

0 

" a 020 


1? 

0 

-* 2 ’ 

0 

a 020 

0 

= 

0 

i 2 

0 

-a 

_ 020 

0 

a 002 


I * 2 

0 

X 2 


(in which the fact was used that = -a 020> as " be no 't e< 3- from 

table l). The determinant of this matrix is 
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A = a 020 ( a 200 a 002 " “OSO 2 ) 


( 9 ) 


and the inverse of this matrix is 


[ M u] s 


Moo 

**01 

Uq 2 

- 

r ' ^02 


n 

a 020 



a 200 a 002 - 

“020^ 


a 200 a 002 _ 

a 020 2 

MlO 

Mil 

M 12 

= 

0 


1 

a 020 

0 


m 20 

M 2 1 

m 22 


a 020 


n 

^OO 



a 200 a 002 “ 

„ 2 
a 020 


a 200 a 002 “ 

a 020 2 


(10) 


Both the determinant and the nonzero elements of occur in 

equation (3 ) . 


Of the 27 third moments only 10 are distinct, and two of those are 
zero as a result of stationarity . The remaining 8 are listed in table 1, 
and expressions are given for them in terms of the triple correlation 
function and double spectrum defined by 


*xx( t :dT 2) = x(t)x(t+T 1 )x(t+T 2 ) 


( 11 ) 


and 


9xx(a>L»“2) = f f e 1 ^ Tl_KU2T ^ 1 l f xx( T l^ T 2) dT l dT 2 ( 12 ) 


These expressions can be derived readily from the definitions of the 
moments , of the correlation function, and of the spectrum by a straight- 
forward extension of the procedure used for the second moments; these 
derivations were obtained by using integrations by parts, differentia- 
tions under the integral sign, and similar elementary operations, and 
by taking advantage of the fact that the process is stationary. Some 
of these moments have been calculated in reference 8 for the purpose 
of est ima t ing the frequency distribution of x(t) alone. 
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Similarly, of the 64 fourth moments the 13 distinct and nonzero ones 
are given in table 1, as are expressions for them in terms of the quad- 
ruple correlation function and triple spectrum 


^XXx( T l> T 2> T 3) = X (t )x (t+T-j) X (t+T 2 ) X (t+Tj) 


(13) 


and 


'Pxxx (01,06,05) s ^3 / M M 6-1 ^ Tl+ ^ T2+a ^ T3 ^ *30* ( t 1 ,t 2 ,t 3) dT l dT 2 dT 3 


(14) 


Expressions for the Expected Number of Maxima 

If the expression for the joint frequency distribution given by 
equation (4) is substituted into equation (2), the following expression 
is obtained for u(xq) : 


n(Xo) 


jAJk, ,] 

1 31 Ml 


where 


- qW 1 jiJB,. , ^ lOq^V™ T ijklmn 


41 


( x o) 


61 


M 


(15) 


I 1 Jk- "( x o) ■/ |* 2 | 

U -00 


a r 


ax^xj^xt . . . 


f G (xo,0,x 2 ^ (2rt)' 5 ^ 2 l/& dx 2 ( 16 ) 


where r is the number of indices i,j, k, .... 

As in the case of the moments, another definition of these functions 
is more convenient for some purposes, namely, 
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3 ™-p ( x o) 


/■° w _i£S2_ 

Sx 0 m a Xl n ax 2 P 


3 /2 

f Cr( x 0 ,0,x 2) (^^ dx 2 


(IT) 


These functions can he identified with those defined by equation (l 6 ) in 
the manner indicated for the moments defined by equations ( 5 ) and ( 6 ), 
respectively. 

These functions are listed in table 2 , in terms of the dimensionless 
variable 


and the parameter 



a 


v/A 

Tg 3 / 2 

x 


(18) 


( 19 ) 


The functions E(z) and E*(z) in this table are defined by 


and 


# 
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where erf and erfc designate the error function and the complementary- 
error function, respectively. The functions not listed in table 2 are 
zero or not required because they are multiplied by a moment which is 
zero. Also included in table 2 for the sake of completeness is the- 
function for Gaussian distributions ^OOO* 

In terms of these functions, equation ( 15 ) can be written as: 


(2j0 5 / 2 ]|a n(xq) » Iqoo - jjf (“300 I 300 t 3a 201 I 2Cl + 3 “l20 I 120 + ^UOehoS + a 003 I 003}j + 

j^T [(^00 " 3a 200 2 ) ■'loo + ^(“soi + 3a 200 a 02o) *301 + f( a 220 ~ a 200 a t)2o) I 220 + 

6 (“202 " ^ OSO 2 ■ a 200 a 002 ) I 202 + 12 ^“121 + a 020 2 ) I 121 + ^(“103 + 3 a 020 a 002 ) I 103 + 
(°t>40 " 5a 020 2 ) 1 0l0 + 6 (“022 “ a 020 a tXI2) I 022 + ( a 004 3a 002 2 ) -tol] + 

Ifb 00 2l 600 + 6a 300 a 201 I 501 + 6 a 300 ct 120 1 l 20 + ( 6a jo0 a a02 + ^Ol 2 ) I 402 + 
l 8 a 201 a 120 1 321 + ( 2a 300 a 003 + l8a 102 a 20l) I 303 + 9a 12Q £l 240 + ( lSa 120“l02 + 

36a lll 2 ) I 222 + ( 6a 201 a 003 + 9ot 102 2 ) I 20l + ^^Ohur 1 !!! + ( 6a O03 a 120 + 

3^ a 012 a lll) X 123 + 6a 102 a 003 I 105 + a 030 2l 060 + 6 ct 030 a 012 iol 2 + 

9a 012 2l 024 + a 003 2 l 006|] - • • • (2°) 

A corresponding expression for the expected maxima per unit time 
exceeding the level x Q can be obtained by integrating this expression 

over xq, inasmuch as 


N 


p 00 

( x °) =J n ( x o') dx o' 

x 0 


(see eq. (l)). The resulting expression for N^Xq) has a form identical 
to equation ( 20 ) but with all functions Xmnp^x^ replaced by func- 
tions Jmnp (xQ) defined by 


p 00 

J mnp( x o) = J -Wrp ') <ix 0 ' 


( 21 ) 


These functions are also given in table 2. Finally, an expression for the 
total number of maxima per unit time can be obtained by replacing Inm.p( x o) 
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in equation (20) by Jmnp ( -“) • For tbe sake of convenience the values 
of Jmnp ( — eo ) ■ are also listed in table 2. 


Input -Output Relations 

If x(t) is the output of a linear system subjected to an input 
F(t ) , the double and higher order correlation functions and spectra 
of x(t), which are required to obtain the moments when no other informa- 
tion is available, can be related to the correspondingly defined function 
for F(t) . In this process either the indicial response h(t), that is, 
the response to an Impulsive input, or the transfer function H(o>), that 
is, the complex amplitude response to steady-state sinusoidal oscilla- 
tions of unit amplitude, must be known. These two characteristics of 
the system are related by 


H(m) = / e“ iojfc h(t)dt (22) 

J 0 


The double correlation functions are related by the expression 


where 


t x ( T ) = [ t F (r-cr)t h (ff)<lcr 

'J -00 





h(t)h(t+ 1 ai )dt 


( 23 ) 


(210 


and the corresponding spectra are related more simply by 

cp x (o>) = | H (co) | 2 q>p(ai) (25) 

where the vertical bars on H(cd) designate the absolute value. 

Similar relations can readily be derived for the higher order 
correlation functions and spectra, the expressions for the spectra being 


?xx K^) = H (^) H (“2) H * (°^ + “2) 9 ff (®L*“2) 


(26) 
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cpxxx^^?^) = H ^ H H ^ H* ^+0^+0^) cpjpF (0^,0^, a>^ (27) 


and so on, where the asterisk designates the complex conjugate. 


DISCUSSION 


Tn the preceding presentation, attention has been confined to maxima. 
However) the results can readily be modified to apply to minima as well, 
because the expected number of minima per unit time^in a band of unit 
width centered on Xq is given by 


n 


p 00 

'( x °) =J o x 2 f(x 0 ,0,x 2 )dx 2 


(28) 


A comparison of this equation with equation (l), which may be written 

as - — - 


n 



x 2 f(Xo,°,-x 2 )dx 2 


(29) 


indicates that the expected number of minima can be obtained for the 
expected number of- maxima by changing the sign of x 2 in the frequency- 

distribution function. 

Consequently, as a result of the definition of 3^, (xq) and 
Jmnp ( Xq) > the expected number of minima can be obtained by making the 
following two changes in the results presented in this paper for maxima: 

(1) Replace M Q2 by (-M Q2 ) wherever it occurs in table 2. 

(2) Multiply those of the functions I n mp( x o) and J mnp( x b) for 
■which n is odd by (-1). 

S imi larly, if the number , of maxima per •unit time below Xq is 
desired, this number can be obtained by subtracting N(xq) from N(-oo) ; 

and in view of the statements made in the preceding paragraph, the same 
procedure can be used for the number of minima below Xq if N^Xq^ is 
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calculated for minima by changing the sign of Mq£ and of the specified 

functions. The total number N(-oo) of maxima is also the number of 
minima, because for each maximum there must be a minimum . 

As may be gathered from the procedure outline herein, when the 
joint frequency-distribution function is not Gaussian, the computational 
effort involved in obtaining the desired expected number of maxima or 
minima soon becomes quite extensive, particularly in view of the fact 
that usually multiple correlation functions or spectra have to be cal- 
culated first in order to obtain the required moments. Although equip- 
ment exists to measure such functions directly or through analog-computing 
devices (see ref. 9 , for instance), this equipment has not reached the 
perfection of the equipment used for the measurement of ordinary correla- 
tions and power spectra. 

The numerical calculations of these functions from time histories 
also pose difficulties beyond those resulting from the greater number 
of variables involved. The source of these difficulties may be described 
by expressing the power spectrum corresponding to the nth correlation 
function associated with a given process x(t) in terms of the Fourier 
transform of the process 


as 



e -icjfc x (t)dt 


(30) 






ltm 
T -» 00 


a(ta I ;T) a (g^;' r ) 


a K-lJ T ) a * K+<V'- * * +< %-l ;T ) 
21 


(51) 


where the asterisk designates the complex conjugate. Hence, in the 
ordinary spectrum (n = 2) only the absolute magnitude of a(co;T) occurs, 
whereas in all higher spectra the phase enters as well. (See also 
ref. 8.) Also, inasmuch as the lim a(o>;T) does not exist, the exis- 

T 00 

tence of the spectra stems from the fact that the product of the func- 
tions a(o;T) is divided by T before the limit is taken. Consequently, 
the higher the spectrum the more the effect of the single T in the 
nominator is "diluted." Consequently, in a practical case, for a given 
reliability, a greater portion of a random process must be analyzed if 
higher-order correlation functions or spectra. are to be obtained than 
if only the ordinary functions are of interest. (See also ref. 10.) 
However, if the deviation from a Gaussian distribution is not too large, 
the higher order moments are likely to be small compared with the second 
moments; therefore, a lower degree of reliability in their determination 
may be acceptable.. 
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In view of these difficulties the procedure outlined .herein is 
likely to find application primarily in problems -where the expected 
number of maxima of an output have to be estimated from knowledge con- 
cerning the input, because then even the statistical! characteristics of 
the process are not known. For a given process the. procedure outlined 
herein is likely to become of advantage only once the equipment for 
measuring higher order correlation functions and spectra is perfected, 
because if the process itself is known. (in the form of a time history) 
the expected number of . maxima can be established by direct count more 
readily than indirectly througi a calculation of the joint frequency 
distribution or moments; and if the frequency function is known., a 
direct calculation of the extrema from equation ( 2 ) or ( 28 ) is likely 
to be more convenient_than an indirect calculation based on moments 
calculated from this function by using equation ( 6 b)". 


CONCLUDING REMARKS 


A method has been outlined for calculating the expected number of 
maxima or minima of a random process with non-Gaussian frequency distri- 
bution from the statistical moments of the process n.nfl its first two 
derivatives. This method is based on an estimate-ofLthe joint frequency 
function of the process and its first two derivatives by means of a 
generalized form of Edgeworth's series; the procedure consists essen- 
tially in applying a correction to the results for a Gaussian process. 
The functions required in this procedure have been calculated for the 
first two correction terms; therefore the effects of skewness and 
kurtosis can be calculated, provided the third and fourth moments of 
the process and its first two derivatives are known. 

Expressions have been given for these moments in terms of multiple 
correlation functions and multispectra, and the relations between these 
functions for the random output of a linear system and those for the 
random input have been indicated . ... 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 19, 1957* 
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TABLE 1.- EGRESSIONS FOR THE MOMENTS IN TERMS OF CORRELATION FUNCTIONS AND SPECTRA 


Moment e 

Expressions in terms of + x (r) , 't r s = c( T 1 > T 2 ^ 5131 

*mcx( t i> t 2' t 3) 

W 

Expressions in terms of <p x (co) , <P xx ( a: l> t “^> end 

*xxx(‘V <a 2’ a> 5) 


“200 

*x 

/"ox** 

Second 

“101 

" 2 “ ~°020 
or 

"“020 


“020 

af** 

&T 2 

pOO 

/ Cpxfl^dai 

^0 


“002 

sSjt 

r- k 

/ qi^ardco 
J 0 


“300 

*xx 

'*» ** ** 


“201 

3ti 2 

- < Pxx?>L 2 'toL ang 


“120 

“ - i “aoi 

- 1 “201 

Third 

“in 

35+„ 

- f J J qwm3“g aa’L ang 


“102 

^♦xx 

Sr^arg 2 

/" f 2<n z Z/iai i ang 

-00 J —CO 


“030 

a^xx 

- 2 aTi&T2 2--^H 

_at lll 


“012 

_ a5 +xx 

- 2 71,/ 93a®i J a>2 2 ao^ dog 


“003 

./*= , *= \ 

- | f [ 9«(aiW'KDL 4 <e2 2 )aa> L dag 

O _«0 


“The functions * x (t), *xx(?i’ t 2)* eni3i +xxx( T i ’ T 2’ T 3) aDi their derivatives are to he evaluated at 
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to-ht.t* i._ E mt E33I0HS TOE IHE MCMEHTS XH THWS OF CORRELATION FUNCTIONS AND SEBCTSA - Concluded 


Moments 

Expressions in terms of tx( T )> ♦xx( T i*’i‘ 2 ), and 

* XXX ( t 1> t 2' t 3) 

(a) 

Expressions in terms of <Px(cd), 9xx( 0J l> a3 2)3 aai 
c Pxxx( <D l7 aj 2J a ^) 


“400 

^XXX 

gj" J J Vyacx *“1 ^ ^ 


“301 

S 2 *^ 
Sri 2 . 

-gj J J cp^j^tOLSdoj. diD2 da? 


“220 

S^JOOC 1 

^ 2 -'?“ 501 

- 7 “301 


g au 

^♦xxx 

- f/JX <P ™^- 2 “2 ^ ^2 ^ 


“202 

aN'aoc 

| [ [ f VjaafrSt**}. *“2 ^ 

«-/ _aa — eo — <0 


*130 

S 5 +xxx 

aT i aT 2 &T 3 

-“on 

Fourth 

“121 

1 3^+XCC 1 .. 

ar^a^ ' 6 -a^T " 9 “202 



“112 

aTl^r^Sr, 

I f f f WYSS ^2 a “5 


“103 

&6 +«x 

^ t 1 2 ^ t 2 2 ^ t 3 2 

-f 1 j" <I) ^o t a5 L 2m 2 2<r 5 S d®L daj, d^ 


“oto 

- 5 a 7^T } ' " 5 “ 121 

■5 a 121 


“022 

f 3 + “« [ „ S \ 
Vl 2 ^ 2 * 1 ^ 2 

J/J1L ^ to 2 ^ 


“013 

a7 

^ t 1^^ t 2 2 ^ t 3 2 

/— to «J /•• 3 2 2 

^tJJJ “2 “2 ^ ^ 


“ocfc 

3 / — j-jga — + — s8 *«« \ 

\ari ! ‘&r 2 2 aT 3 2 &ri 3 aT 2 5 ar 3 2 y 

| J J J (j°i ^ ^nig dm, 


»The functions y x (r), aid * Jncc (T 1 ,T 2 , T 3 ) and their derivatives ere to be evaluated at 


t - T ! ” T 2 - T 3 


- 0 . 
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